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Report  Abstract 


The  two-year  subject  program,  conducted  through  tight  coupling  between  experiment,  theory 
and  computation,  and  reported  in  high  impact  journal  articles,  has  focused  on  the  dynamics  of 
turbulent  flames  mainly  for  hydrocarbon  fuels  in  environments  simulating  various  operational 
aspects  of  aero-engines.  The  thrust  for  this  program  constitutes  of  three  major  areas  of  turbulent 
combustion:  (1)  Flame  surface  statistics,  (2)  Flame-turbulence  interaction,  and  (3)  Effects  of 
turbulence  on  the  growth  of  ignition  kernel. 

Using  high-speed  Mie- scattering  we  investigated  the  local  flame  surface  statistics  of  constant- 
pressure  turbulent  expanding  flames.  First,  the  statistics  of  the  local  length  ratio  is 
experimentally  determined  from  high-speed  planar  Mie  scattering  images  of  spherically 
expanding  flames,  with  the  length  ratio  on  the  measurement  plane,  at  predefined  equiangular 
sectors,  defined  as  the  ratio  of  the  actual  flame  length  to  the  length  of  a  circular-arc  of  radius 
equal  to  the  average  flame  radius.  Assuming  isotropic  distribution  of  such  flame  segments  we 
then  convolute  suitable  forms  of  the  length-ratio  probability  distribution  functions  (pdfs )  to  arrive 
at  the  corresponding  area-ratio  pdfs.  It  is  found  that  both  the  length  ratio  and  area  ratio  pdfs  are 
near  log-normally  distributed  and  shows  self-similar  behavior  with  increasing  radius. 

To  quantify  the  effect  of  the  flame  on  the  turbulence  the  mean  and  the  rms  of  the  flow  field  are 
measured  over  a  range  of  conditions  both  in  the  presence  and  absence  of  the  flame  using  High 
Speed  PIV.  The  distributions  of  stretch  rate  contributions  from  different  terms  such  as  tangential 
straining,  normal  straining  and  curvature  are  also  obtained.  It  is  found  that  the  normal  straining 
displays  non-Gaussian  pdf  tails  whereas  the  tangential  straining  shows  near  Gaussian  behavior. 
We  have  further  tracked  the  motion  of  the  edge  points  that  reside  and  co-move  with  the  edge  of 
the  flame  kernel  during  its  evolution  in  time,  and  found  that  within  the  measurement  conditions, 
on  average  the  persistence  time  scales  of  stretch  due  to  pure  curvature  exceed  those  due  to 
tangential  straining  by  at  least  a  factor  of  two. 

Finally  we  experimentally  demonstrated  that  contrary  to  the  belief  that  ignition  of  a  combustible 
mixture  by  a  high-energy  kernel  is  more  difficult  in  turbulence  than  in  quiescence  because  of  the 
increased  dissipation  rate  of  the  deposited  energy,  it  can  actually  be  facilitated  by  turbulence  for 
mixtures  whose  thermal  diffusivity  sufficiently  exceeds  its  mass  diffusivity.  In  such  cases, 
turbulence  breaks  the  otherwise  single  spherical  flame  of  positive  curvature,  and  hence  positive 
aerodynamics  stretch,  into  a  multitude  of  wrinkled  flamelets  possessing  either  positive  or 
negative  stretch,  such  that  the  intensified  burning  of  the  latter  constitutes  local  sources  to 
facilitate  ignition. 
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1.  Experimental  Setup: 


All  the  experimental  studies  reported  here  used  a  constant-pressure,  fan-stirred  combustion 
chamber  to  investigate  the  propagation  of  a  spherically  expanding  flame  (Fig.  1.1).  Chambers 
based  on  this  design  have  been  used  to  study  various  aspects  of  laminar  and  turbulent  expanding 
flames,  and  are  detailed  in  our  previously  published  articles  [Tse  et  al  (2004),  Kelly  and  Law 
(2009),  Chaudhuri  et  al  (2012)].  Briefly,  the  apparatus  consists  of  two  cylindrical  chambers:  an 
outer  chamber  that  houses  an  inner  chamber  of  a  much  small  volume  and  an  almost  unity  aspect 
ratio  (Inner  dia:  114  mm;  outer  dia:  167  mm;  length:  127  mm).  The  flat  ends  of  the  inner  chamber 
are  sealed  with  thick  quartz  windows  for  optical  access.  The  inner  chamber  is  filled  with  the  test 
fuel/air  mixture  which  is  spark-ignited,  while  the  outer  chamber  is  filled  with  an  inert  of  the 
same  pressure.  The  connection  between  the  chambers  is  controlled  by  four  valves  (sliding  sleeve 
with  multiple  holes)  that  are  opened  at  the  instance  of  the  spark  discharge,  hence  allowing  the 
expanding  flame  to  come  in  contact  with  the  inert  in  the  outer  chamber  and  instantly  quench.  The 
expansion  process,  thus,  is  nearly  isobaric  owing  to  the  small  volume  of  the  inner  chamber 
compared  to  that  of  the  outer  chamber.  Experiments  up  to  the  range  of  30  to  60  atm  have  been 
conducted  based  on  this  dual-chamber  design,  although  the  maximum  pressure  for  the  current 
experiments  was  kept  to  5  atm  as  it  is  sufficient  for  the  present  study. 


Experimental  Setup.  CV:  Check  Valve,  PG:  Pressure  Gauge,  PT:  Pressure 
Transducer,  M:  Fan  Motor,  L:  Cylindrical  Lens,  E:  Electrodes,  W:  Quartz  Window 


Fig.  1.1:  experimental  setup  with  all  the  components. 


The  turbulence  in  the  chamber  is  generated  using  four  equally  spaced  continuously  running  fans 
(69  mm  diameter),  powered  by  motors  enclosed  in  the  outer  chamber.  The  motors  are  controlled 
by  using  a  feedback  mechanism,  which  regulates  the  power  of  the  motors  to  maintain  constant 
speed. 
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To  detect  the  flame  edge  during  the  combustion  process  we  seed  the  air  using  DEHS  droplets 
(nominal  diameter  ~  1-2  jum)  generated  by  a  high-pressure  nebulizer.  A  high-speed  Nd-YLF  laser 
(527  nm)  synchronized  with  a  high-speed  camera  (Phantom  v7.3)  is  used  for  Mie  scattering. 
Using  beam  expanding  optics  the  2  mm  diameter  laser  beam  is  expanded  to  a  750  jum  thick  laser 
sheet  inside  the  inner  chamber  as  shown  in  the  figure.  The  camera  placed  at  90°  from  the  laser 
plane  recorded  the  event  at  8000  fps  with  a  resolution  of  608x600  to  achieve  a  spatial  resolution 
of  -50  /mi/ pixel.  Before  ignition  the  entire  chamber  remains  filled  with  droplets  showing  a 
uniform  distribution  of  seed  density.  Following  ignition,  as  the  flame  propagates  the  seeding 
droplets  inside  the  flame  vaporize,  hence  limiting  the  high  seeding  density  only  outside  the  flame 
(Fig.  2bi).  Post  processing  the  high-speed  images,  the  instantaneous  2D  flame  (on  the  laser 
plane)  edge  (the  edge  of  the  seeding  concentration)  was  detected  using  Canny  edge  detection 
technique  through  Matlab.  During  image  processing  we  mask  the  electrodes  to  avoid  their 
influence  of  the  statistics. 

2.  Flame  Surface  Statistics: 

Turbulent  flame  speed  is  a  major  topic  of  sustained  interest  in  turbulence  and  combustion 
research.  Starting  with  Damkohlcr's  pioneering  work  [Damkohler  1940],  it  is  generally  believed 
and  shown  by  many  studies  that,  at  least  for  equidiffusive  mixtures,  the  turbulent  flame  speed  is 
largely  controlled  by  the  total  flame  surface  area,  with  the  normalized  turbulent  flame  speed 
assumed  proportional  to  the  normalized  flame  surface  area.  Since  the  turbulent  flame  speed  is  a 
global  statistics,  it  is  of  interest  to  consider  how  the  distribution  of  quantities  related  to  the  local 
flame  surface,  such  as  the  flame  surface  area  ratio,  build  up  this  global  measure.  The  turbulent 
flame  surface  is  typically  highly  wrinkled  and  folded  over  a  multitude  of  scales,  starting  from  the 
hydrodynamic  length  scale  to  the  smallest  scales  controlled  by  the  Kolmogorov  length  scale  and 
flame  properties  such  as  the  Markstein  length  and  flame  thickness.  Therefore  it  is  useful  if  the 
information  contained  in  this  complex  geometry  can  be  projected  onto  a  simple  regular  geometry 
such  as  a  plane  or  a  smoothly  curved  surface.  A  straightforward  way  to  arrive  at  this  simple 
geometry  is  to  just  consider  a  filtered/average  location  of  the  flame  positions  to  represent  a 
smooth  surface.  The  information  contained  in  the  wrinkled  surface  if  projected  on  this  simple 
geometry  will  also  make  it  amenable  for  well-known  processing  techniques  such  as  spectral, 
wavelet  or  multifractal  analysis.  Not  only  such  an  approach  would  be  useful  for  analysis  but  it 
might  be  also  of  interest  in  reduced-order  modeling,  such  as  Farge  Eddy  Simulation  or  Reynolds 
Averaged  Navier-Stokes  Equation,  in  which  very  fine  scale  flame  surface  wrinkling  and  the 
information  contained  therein  cannot  be  directly  simulated  due  to  limitation  of  the  resolution.  By 
this  approach  the  wrinkled  surface  information  randomly  distributed  in  three  dimensions 
(volume)  is  now  described  using  two  dimensions  (a  regular  surface),  essentially  resulting  in  a 
reduction  of  the  dimensionality  of  the  problem.  The  knowledge  of  the  local  area  ratio  is  also 
pivotal  in  obtaining  the  Probability  Distribution  Functions  (pdfs)  of  any  quantity  conditioned  on 
that  surface.  Motivated  by  its  importance  for  both  fundamental  understanding  and  reduced-order 
simulation  of  turbulent  flames,  this  work  focuses  on  the  local  statistics  of  flame  surface  area. 
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Ensemble  averaging  of  flame  edges  relative  to  their  centers 


(a) 


Fig.  2. 1  (a)  This  figure  shows  instantaneous  realizations  of  the  flame  edges  (blue  lines)  at  the  same 
<R>  for  all  the  conditions  and  runs.  The  thick  red  line  shows  a  curve  obtained  by  ensemble 
averaging  of  instantaneous  flame  edges  at  each  angular  sector,  (b)  This  image  demonstrates  the 
procedure  of  extracting  the  length  ratio  statistics  from  Mie  scattering  images,  (i)  A  typical  Mie 
scattering  image.  Red  line  shows  the  flame  edge,  green  line  shows  the  circle  with  mean  radius  <R>C. 
(ii)  Extracted  flame  edge  broken  down  to  angular  sectors,  (iii)  In  an  angular  sector  the  important 
measurements,  dl0,  dlc  and  dlr  It  is  important  to  note  that  in  this  sector  the  flame  edge  overcuts  the 
mean  circle  ( i.e .,  the  flame  edge  is  at  the  higher  radial  location).  However,  in  some  angular  sectors  it 
can  undercut  the  mean  circle.  Thus,  the  quantity  (LR) j  can  locally  have  values  less  than  unity. 


Using  the  G-equation  and  starting  with  the  flame  surface  area  ratio,  Kerstein  et  al.  (1988) 
obtained  an  expression  for  the  turbulent  flame  speed  of  a  statistically  planar  flame,  propagating 
in  the  i  direction  in  the  mean.  The  projected  area  of  such  a  propagating  flame  surface,  i.e.  the  y-z 
plane,  was  divided  in  equal  infinitesimal  areas  d\  =  dydz .  Consequently,  if  the  local  flame 

surface  area  contained  within  a  cylinder  of  cross  sectional  area  d\  with  the  axis  parallel  to  i  is 
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given  by  JS(G0),  where  i  is  the  unit  normal  in  ^-direction  with  the  flame  being  represented  by  a 
level  set  G  =  G0,  then  the  turbulent  flame  speed  can  be  expressed  as: 


ST  \  /dS(G„)\ 

SL  \  \  dA,  j 

It  was  recognized  that  due  to  folding  of  the  flame  surface  G  =  G0 ,  multiple  flame  segments  might 
be  encountered  within  a  given  cylinder.  Each  of  these  segments  referenced  by  index  j  was 
considered  locally  planar  with  local  surface  normal  represented  by  n  .  Then  the  area  ratio  for  the 

j th  segment  is  given  by  dS  (G0  )  /  dA{)  =  1  /  |i-ny  |  and  the  flame  area  contained  within  the  cylinder 
is  given  by  summing  over  j  : 

J5(G0)/cM0=^l/|i«n.|.  (2.2) 

j 

Correspondingly,  for  an  expanding  flame  which  is  nearly  spherical  (from  experiments  reported 
in  the  following)  on  ensemble  averaging,  as  shown  in  Fig.  2.1a,  a  sphere  with  radius  equal  to  the 
average  radius  of  the  flame  can  be  considered  to  be  the  base  or  projected  surface.  This  projected 

surface  can  be  decomposed  into  infinitesimal  area  elements  dAs  =  {R)*  sin  6d0d(f> .  where  ( R ) 

is  the  average  radius,  the  polar  angle  (j)  and  the  azimuthal  angle  9.  To  assess  the  ratio  of  the  actual 
flame  surface  with  the  projected  area,  we  consider  imaginary  tetrahedral  cones  whose  apex  is  the 
center  of  the  sphere  and  has  cross-sectional  area  dAs  where  it  intersects  the  projected  surface  (as 
shown  in  Fig.  lb).  The  wrinkled  flame  area  contained  within  a  cone  can  be  represented  by 
dSs  (G0) .  In  our  experiments,  measurements  are  performed  using  2D  Mie-scattering  technique  at 

the  laser  plane,  hence  the  above  discussion  should  be  modified  for  a  planar  case  where  the 
projected  contour  is  a  circle  and  the  infinitesimal  projected  length  element  is  given  by 
dlQ  =  (R)^  d(f> .  In  the  above  the  subscripts  5  and  c  represent  spherical  (3D)  and  circular  (2D) 

projections  respectively.  The  corresponding  length  ratio  contained  within  each  of  the  sectors  is 
given  by 

dl(G„)ldl0=Y.(il\v,nJ\)(r1/(R}c).  (2.3) 

i 

Here  r;  is  the  unit  normal  in  the  radial  direction  and  the  second  factor  arises  due  to  change  in 
length  of  the  circular  segment  between  the  local  radius  r  and  average  radius  (/?)  .  Summation 
over  j  is  still  important  as  the  flame  edge  determined  can  be  folded,  resulting  in  multiple  flame 
segments  crossing  within  a  given  sector.  Schematic  of  the  determination  of  such  a  length  ratio  is 
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shown  in  Fig.  lb.  Now  the  flame  edge  that  we  obtain  can  be  represented  by  the  level  set  G  =  G0 , 
with 


n  =  VGg_g 

]  O-lr0 


/  VG 


lG=Gn 


Hence  Eq.  (2.2)  can  be  rewritten  as: 

dl(G„)l dl„  =  Z(| VG|0=0i  / |r.VG|c=0< )(r.  /  (*), ) 

j 

which  in  r  -  f  coordinates  is  given  by, 


(2.4) 


(2.5) 


dl(G0)/dl0=XJ  1- 


^isg^2  /r 


8G" 
\dr  y 


rj 


(2.6) 


z.e.  in  this  2D  realization  any  surface  fluctuation  in  the  out-of-plane  direction  (here  6)  is 
neglected.  In  particular  the  actual  area  ratio  in  the  r-(f>-6  spherical  coordinates  is  given  by: 


dS(G0)/dA0  =J\| 
J  V 


1+ 


A1  SGV 


r  d(j) 


' 8G 
v  dr 


1  dG 

rsin^  86 


' 8G 
ydr  y 


rj 


WRlsJ 


(2.7) 


We  will  discuss  shortly  that  using  statistical  methods  and  some  assumptions  about  the 
corresponding  distributions  one  can  arrive  at  the  3D  area  ratio  statistics  from  the  measured  2D 
length  ratio  statistics.  While  extrapolation  from  2D  statistics  to  3D  statistics  is  not 
straightforward,  it  can  be  certainly  expected  that: 

ST  , /dS(Gc)\  /dl(G0)\ 

St  \  d\  I  \  dl„  I 

Recognizing  this,  we  propose  a  systematic  method  for  converting  the  2D  length  ratio  statistics  to 
the  3D  area  ratio  statistics.  In  doing  this  we  assume:  i)  the  statistics  of  the  surface  are  invariant 
upon  rotation;  and  ii)  the  in-plane  and  out-of-plane  statistics  are  mutually  independent. 
Consequently,  we  can  utilize  the  pdf  of  the  second  term  in  the  RHS  of  Eq.  (2.7)  to  generate  the 
pdf  of  the  corresponding  third  term.  Then  we  need  to  generate  a  pdf  of  the  sum  of  these  two 
terms  which  can  be  done  by  convoluting  the  two  pdfs  as  shown  below. 

Let  X  and  Y  be  two  discrete  and  independent  random  variables  with  Z  =  X  +  Y  as  the  sum  of 
the  two  random  variables.  Then  the  events  X  =  a  and  Z  —c  can  only  occur  simultaneously  if 
Y  —  c  —  a  .  Since  X  and  Y  are  independent,  the  probability  of  Z  -  c  is  thus  given  by  convolution 
of  the  probabilities  X  -  a  and  Y  =  c-a , 
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(2.9) 


oo 

P(Z  =  c)=  YJp{x=a)-p{Y=c-a) 

a=-o o 


Through  edge  detection  technique  used  on  2D  Mie-scattering  images,  the  instantaneous  flame 
edge  can  be  extracted  as  shown  in  Fig. 2.  lb  (ii).  It  is  to  be  noted  that  the  flame  edge  at  any 
instance  is  representation  of  the  level  set  function  G  =  Go.  The  figure  also  shows  the  circle 
(centered  at  the  centroid  of  the  area  under  flame  edge)  with  a  mean  radius  (R)  ■  The  flame  edge 

has  been  divided  into  small  angular  sectors  of  l°each.  Through  image  processing  technique  we 
calculate  or  evaluate  the  ratio  of  the  flame  edge  length  to  the  arc  length  of  a  mean  circle 
( dl(Go)/dlo )  as  shown  in  Fig.  2.1.b  (iii).  Due  to  stretching,  wrinkling  and  folding  under  turbulence 
it  is  possible  that  the  flame  edge  can  have  multiple  cuts  or  segments  in  one  sector.  Although  the 
number  of  segments  in  each  sector  can  be  different,  each  sector  will  have  at  least  one  flame 
segment. 

Now,  following  Eq.  (2.6),  for  the  /th  flame  segment  at  the  zth  sector  we  can  write, 


[dl{GQ)/dlQ\,j 


r\ 5GV 

vr  d</>, 


' 8G ' 
V  Sr  , 


(2.10). 


We  next  define  a  parameter  Uj  .  = 


value  of  i  and  j  using 


vA,; 


1  dG 
r  d(j) 


f  dG}2 

V  Sr  j 


which  can  be  evaluated  for  any 


f  \ 

riJ 

4 

r\  dG  ' 

7i 

' dG ' 

2 

r 

(R) 

\\  G ) 

Kr  Sf  J 

r 

V  Sr  J 

V 

R 


[dl(G0)/dlQlj 


'  ^  v 
VV^y 


(2.11) 


•th 


From  the  statistics  of  Ui ; ,  it  is  possible  to  obtain  a  probability  distribution  function  for  the  j 
segment  i.e.  pdf  \jJ  j  J . 

Now  if  we  assume  that  the  statistics  of  the  surface  are  invariant  upon  rotation,  i.e.  the  surface  is 
isotropic,  we  can  say  experiments  on  a  plane  rotated  by  90°  [r-0  plane],  the  statistics  of  the 


quantity 


1  dG 


' dG ' 
v  Sr  y 


would  have  been  similar  to  the  statistics  of 


\r  SO  j 

obtained  for  the  current  experimental  plane  (r-tp  plane)  or 


r\  SGX  /(dG x 

r  dtp  J  /  [dr  y 


Pdf 


'  r.  ' 


1  5G 
r  SO 


dG 

dr 


=  Pdf 


'  r  A 


rl  dCf'2 

r  dp 


dG 

dr 


=  pdf[uj]  (2-12). 


8 


Now  as  we  know  the  statistics  that  sin^  is  based  on,  we  can  also  obtain  pdf\j']  ,  where 


ru 


'  r..  ^ 

vA; 


1  dG 
r  sirup  36 


Following  the  definition  of  U,  we  define  a  second  parameter 


X.  = 

i>J 


A 


R) 


UJ 


1  + 


r\ dG^2 

r  d(j> 


(2.13a). 


(2.13b) 


and  calculate  the  probability  distribution  function  for  the  jth  segment,  i.e,  pdf  \^X  i  J . 

Using  our  second  assumption  that  the  statistics  on  the  r-rp  plane  and  r-0  plane  are  mutually 
independent,  we  can  use  the  convolution  between  pdf\ji^  and  pdf  Y j  J  to  arrive  at 

.&] 

for  all  possible  values  of  j,  we  arrive  at 


pdf  j  +  Yj  ] .  Now  we  define  Z;  =  AX  +  Y] ,  and  perform  successive  convolution  on  pdf 


Pdf 


dS(G„)/d\  =  Y.Jz~ 


(2.14) 


In  the  following,  we  first  briefly  discuss  the  length  ratio  pdfs  and  then  the  area  ratio  pdfs  in  more 
detail,  for  unity  Lewis  number  ( Le ),  CFL-air,  ^  =  0.90  flames  at  two  levels  of  turbulence 
intensity  s  =  1.43  and  2.85  m/s  with  fan  speed  =  2000  and  4000  rpm,  respectively)  and  at 
three  different  pressures  (p  =  1,2  and  5  atm).  The  Markstein  numbers  {Mk)  obtained  from 
laminar  expanding  flame  experiments  for  mixture  conditions  corresponding  to  the  present 
turbulent  experiments  are  nearly  constant,  slightly  decreasing  with  pressure  from  1.18  to  0.73. 
We  also  include  a  Le  >  1  case:  C4Hi0-air,  <f  =  0.8  (Wnns  =  1.43  m/s  (fan  speed  =  2000  rpm),  p  =  5 
atm,  Mk  =  2.75)  to  qualitatively  assess  the  importance  of  Mk  on  the  flame  surface  area  ratio 
statistics. 

We  next  study  the  flame  length  ratio  obtained  directly  from  the  experiments,  without  any 
assumption.  As  explained  earlier  (Eq.  2.8)  the  length  ratio,  ( LR=dl(Go)/dlo )  is  not  equivalent  to 
the  area  ratio,  ( AR=dl(Go)/dlo ),  which  is  directly  related  to  the  flame  speed.  Nevertheless,  the 
mean  of  this  scalar  quantity  should  follow  similar  characteristics  as  the  flame  burning  flux  which 
increases  with  pressure,  turbulence  intensity,  and  in  general  with  the  turbulence  Reynolds 
number.  Moreover,  with  the  current  experimental  procedure  being  2D  in  nature,  the  length  ratio 
is  the  most  fundamental  parameter  of  interest,  denoting  the  projection  of  the  2D  flame  edge  on 
the  mean  circle.  Since  turbulence  wrinkles  the  flame  to  destroy  its  smooth  spherical  shape  that 


9 


one  can  observe  in  laminar  experiments,  the  length  ratio  is  thus  a  direct  measure  of  the  extent  of 
flame  wrinkling. 
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Fig.  2.2:  Evolution  of  length  ratio  pdfs  with  the  flame  radius  <R>  (mm)  for  CH4  flame:  a) 
latm,  2000 r pm,  b)  1  atm,  4000  rpm  and  c)  Evolution  of  <LR>  with  the  flame  radius  <R> 
(mm)  for  different  CHd  flame  at  different  pressure  and  rpm 


Figure  2.2  shows  the  probability  density  functions  (pdfs)  for  two  different  conditions  at  different 
mean  radii.  The  first  observation  that  emerges  from  these  pdfs  is  that  they  display  a  somewhat 
distorted  log-normal  nature.  Physically  such  a  distribution  signifies  that  the  length  ratio  can 
attain  extreme  values  with  wide  separations,  which  is  an  inherent  nature  of  a  dissipative  quantity 
under  turbulence,  due  to  its  multiplicative  randomness.  Secondly,  it  is  also  evident  that,  despite 
its  similarity  in  the  shape,  the  distribution  at  any  condition  shifts  to  the  right  with  the  radius.  This 
shift  can  be  captured  more  effectively  with  the  mean  LR  (<LR>)  plotted  in  Fig.  2.2e,  which 
shows  a  clear  increase  with  radius.  This  increase  signifies  that  the  flame  becomes  more  and  more 
strained  and  wrinkled  as  is  expands,  which  in  turn  ensures  an  increase  in  the  turbulent  flame 
speed.  Higher  turbulence  intensities  increase  the  small  structures  in  the  flame  edge,  making  the 
flame  more  wrinkled.  This  is  reflected  as  the  increase  in  <LR>  when  we  compare  the  CH4 
flames  under  the  same  pressure  but  two  different  turbulence  intensities  (Fig.  2.2e).  Although  the 
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turbulence  promotes  wrinkling  and  stretching  of  the  flame  edge,  they  are  also  being  inhibited  by 
Markstein  diffusion  which  prevents  the  evolution  of  small-scale  structures.  The  effect  of 
Markstein  diffusion  and  its  influence  on  turbulent  flame  wrinkling  is  not  very  clear  from  the 
<LR>  values  of  the  CH4  and  C4H10  flames  at  5  atm,  2000  rpm,  where  the  latter  has  larger 
Markstein  numbers.  The  plot  (Fig.  2.2e)  shows  that  both  flames  have  very  similar  length  ratios 
making  it  difficult  to  distinguish  the  effect  of  the  Markstein  number.  We  recognize  that  at  higher 
pressures  (such  as  5  atm )  the  smallest  scales  of  flame  surface  wrinkling  may  be  near  or  below  the 
limit  of  imaging  resolution,  barring  proper  reconstruction  of  the  local  length  ratio.  Higher 
resolution  imaging  is  necessary  to  clarify  such  differences  at  these  pressures.  We  will  however 
also  show  later  that  this  influence  of  the  Markstein  number  becomes  clearer  when  we  compare 
the  area  ratio  statistics. 

We  also  have  computed  the  area  ratio  quantities,  as  AR  directly  correlate  with  the  turbulent 
flame  speed.  Figure  2.3  compares  the  Area  Ratio  pdfs  for  two  different  conditions.  Although  the 
shape  of  the  pdfs  for  LR  and  AR  are  similar,  the  mean  of  these  quantities  are  different.  As 
explained  in  Eq.  (2.1),  the  mean  <AR>  is  representative  of  the  flame  speed,  which  is  shown  in 
Fig.  2.3c.  Clearly,  the  values  of  <AR>  are  quite  different  from  those  of  <LR>  presented  in  Fig. 
2.2c.  This  therefore  demonstrates  that  the  2D  measurement  of  the  length  ratio  cannot  readily 
represent  the  normalized  turbulent  flame  speed  before  applying  this  or  other  equivalent 
conversion.  The  values  of  the  area  ratio  presented  in  Fig.  7,  on  the  other  hand  are  obtained  from 
a  3D  statistics,  and  as  such  can  be  used  to  represent  the  normalized  flame  speed.  The 
comparative  study  of  Fig.  2.3c  shows  some  fundamental  characteristics  of  the  flame  surface 
under  turbulence.  As  explained  earlier,  the  flame  surface  undergoes  extensive  stretching, 
wrinkling  and  folding  in  the  presence  of  turbulence.  Due  to  the  multiplicative  nature  of 
randomness,  the  area  ratio  shows  a  near  log-normal  distribution  which  is  a  common  signature  of 
a  dissipative  scalar  in  strong  turbulence.  For  a  turbulent  expanding  flame  the  hydrodynamic 
length  scale  increases  with  flame  size  which  is  an  artifact  of  increased  flame  brush  thickness  or 
the  volume  that  the  wrinkled  flame  occupies.  In  other  words,  as  the  flame  expands,  it  becomes 
more  wrinkled  resulting  in  an  increase  in  <AR>  as  seen  in  Fig.  2.3c.  As  observed  for  the  length 
ratio,  the  area  ratio  also  increases  with  the  turbulence  intensity  due  to  intense  stretching  and 
folding  of  the  flame  [Fig.  7  (CH4  flames:  1  atm  (2000,  4000  rpm )  and  2  atm  (2000,  4000  rpm)]. 
The  effect  of  Markstein  diffusion  on  the  other  hand  becomes  apparent  while  comparing  the  5  atm 
flames  of  CH4  and  C4H10,  in  which  the  former  shows  higher  values  of  <AR>  due  to  more 
wrinkling  of  the  flame  compared  to  C4H10,  which  has  a  higher  Markstein  number. 
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Fig.  2.3:  Evolution  of  normalized  Area  ratio  pdfs  with  the  flame  radius  <R>  (mm)  for  CH4 
flame  (linear  Y  scale  in  inset):  a)  1  atm,  2000  rpm,  b)  1  atm,  4000  rpm,  and  c)  Evolution  of 
<AR>  with  the  flame  radius  <R>  (mm)  for  different  CHj  flame  at  different  pressure  and 


3.  Flame-Turbulence  Interaction 

This  part  of  the  study  focuses  on  the  interaction  between  turbulence  on  flames  and  vice  versa. 
The  fundamental  parameters  that  quantify  the  state  of  turbulence,  if  it  is  nearly  isotropic,  are  the 
turbulence  intensity  unm  and  the  integral  length  scale  Lj.  The  quantities  that  characterize  the  local 

flame  structure  are  the  local  laminar  flame  speed  SL  and  flame  thickness  Sl.  According  to  the 
hydrodynamic  theory  of  premixed  flames  the  stretched  laminar  flame  speed  SL  is  given  by: 

Sl  =  Sl-Sm  K  (3.1) 

where  SL  and  SM  are  the  planar  laminar  flame  speed  and  Markstein  length  respectively  0,  K  is 
the  total  stretch  rate  which  can  be  further  decomposed  as  0: 

stretch  rate  by  curvature:  Kk. 

K=  Slk  -  (v-n)/c  +  Vf-vf  ,  (3.2) 

stretch  rate  by  pure  curvature:  Kc  normal  strain-  K  tangential  strain:  K, 
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in  which  Kt=[dij-njnj^Sij,  and  S:j  =  1  for  /  =  j;  0  for  /;  with  5i;  being  the  strain  rate 

tensor.  The  second  and  third  terms  on  the  RHS  of  Eqn.  (3.2)  together  is  the  stretch  rate 
contribution  by  the  strain  rate,  K s .  Stretch  rates  originating  from  the  tangential  strain  rate  Kt, 
normal  strain  Kn  and  stretch  due  to  pure  curvature  K,  can  be  considered  to  be  parameters  that 
quantify  the  turbulence-flame  interaction  as  well  0.  In  general,  turbulence  results  in  multi-scale 
wrinkling  of  the  flame  surface,  affecting  the  aerodynamics  of  the  flame  which  in  turn  affects  the 
nearby  flow  field  by  heat  release  and  gas  expansion.  Since  the  burned  products  in  the  expanding 
flame  should  be  statistically  stationary,  the  flame  could  induce  a  mean  outward  flow  that  might 
affect  the  turbulence  parameters  like  unns  and  L/  in  the  otherwise  cold  non-reacting  flow.  It 
appears  that  in  all  experimental  efforts  on  turbulent  expanding  flames  only  cold  flow  turbulence 
parameters  like  urms  and  L/  measured  in  the  absence  of  flame  have  been  reported.  It  is  thus 
critical  to  investigate  to  what  degree  these  fundamental  quantities  are  affected,  if  at  all,  by  the 
flame  as  these  modified  parameters  are  what  distort  the  flame  instead  of  the  quantities  measured 
in  its  absence. 


In  view  of  such  a  need,  we  attempted  to  provide  experimental  high-speed  particle  image 
velocimetry  data  that  quantifies  flame-turbulence  interaction  dynamics  well  resolved  in  two 
dimensional  space  and  time,  as  the  flame  evolves  from  a  wrinkled  flame  kernel  to  a  well- 
developed  expanding  turbulent  flame.  The  high-speed  Mie  scattering  which  is  a  necessary 
precursor  to  HS-PIV  allows  us  to  extract  the  flame  edges  from  the  maximum  seed  density 
gradients.  This  enables  tracking  the  evolution  of  the  curvature  statistics.  The  flow  velocity  and 
strain  rate  tensor  components  from  the  HS-PIV  could  be  superposed  on  these  flame  edges, 
allowing  investigation  of  the  evolution  of  stretch  rate  statistics  due  to  both  local  strain  rate  and 
curvature. 

Figures  3.1a  and  3.1b  show  the  mean  radial  velocity  (Ur}  and  the  effective  turbulent  intensity  on 


the  flame  u'eff 


for  all  the  cases  under  study, 


respectively.  We  also  compared  these  results  with  those  obtained  from  the  corresponding  cold 
flow  measurements  in  the  absence  of  flame,  prior  to  ignition.  We  observe  that  in  the  absence  of 
the  flame  the  mean  radial  flow  is  directed  inward  (towards  the  center  of  the  chamber)  indicated 
by  the  negative  values  in  Figure  3.1a.  This  weak  radial  flow  also  increases  with  radius,  closer  to 
the  fans.  However  during  flame  expansion,  the  mean  radial  flow  adjacent  to  the  flame  is  radially 
outward  in  nature  shown  by  the  positive  values  of  (C/r).  As  the  flame  expands  it  pushes  the 
unbumt  gas  in  front  of  it.  This  effect  becomes  stronger  as  the  flame  accelerates  with  increase  in 
radius,  resulting  in  increasing  (Ur )  with  (R\.  However,  we  found  that  there  is  no  discernible 
change  in  the  u'eff  in  the  presence  of  the  flame.  This  shows  the  effect  of  flame  on  turbulence 
intensity  is  rather  weak. 
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Fig.  3.1:  Plot  of  (a)  ensemble  averaged  <Ur>  v,v  <R>  and  (b)  ensemble  averaged  a vs.  <R>  in 


presence  and  in  absence  of  flame 


In  Figure  3.2  we  look  into  the  evolution  of  the  probability  distribution  functions  (pdfs)  for  K,  and 
Kn  defined  in  Eqn.  (3.2).  The  pdfs  are  computed  over  the  ensemble  of  all  the  runs  for  a  given 
condition,  conditioned  at  particular  values  of  (i?) .  It  is  clearly  seen  that  both  the  Kt  and  Kn  show 

self-similar  profiles  during  their  evolution  with  However,  when  plotted  on  a  log  scale  the 

nature  of  the  pdf  tails  of  K,  and  K„  appear  different.  The  pdf  for  K,  shows  parabolic  tails 
suggesting  Gaussian  distribution,  whereas  Kn  shows  non-Gaussian  stretched  tails  implying  finite 
probability  of  very  large  Kn  values  in  comparison  to  its  near  Gaussian  counterpart  K,.  Clearly 
with  increase  in  unns  both  the  K,  and  Kn  pdf  tails  are  flared.  As  the  turbulence  Reynolds  number 
Rej  increases  (with  turbulence  intensity),  the  Kolmogorov  length  scale  //  decreases,  resulting  in 
the  emergence  of  smaller  scale  eddies  which  are  also  highly  strained  and  intermittent.  These 
impart  larger  Kt  and  Kn  observed  from  the  flaring  of  the  tails.  It  is  also  important  to  see  if  the 
increased  unns  or  Rej  has  similar  impact  on  the  mean  stretch  rates  in  addition  to  that  on  the  higher 
moments. 
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Fig.  3.2  :  (a)  PDFs  ofK,  and  (b)  PDFs  ofK„  conditioned  at  specific  values  <R>,  over  six  different 


conditions  of  fuel-air  mixtures,  unns  and  pressure:  (i)  CH4:  latm  and  2000rpm,  (ii)  CHp  l  atm  and 
4000 r pm 


Recently  Clavin  and  GranOtero  0  analytically  showed  the  existence  of  two  unique  Markstein 
lengths,  one  due  to  strain  rate  and  the  other  due  to  curvature.  Following  the  same  notion,  if  one 
can  isolate  the  persistence  time  scales  of  these  two  Markstein  lengths,  it  is  possible  to  distinguish 
between  their  relative  importance  and  to  assess  the  relative  relevance  of  local  strain  and  local 
curvature  for  dissipating  turbulent  flame  surface  fluctuations.  From  a  modelling  point  of  view 
this  is  of  considerable  interest  especially  in  the  context  of  G-equation  used  for  the  description  of 
a  premixed  flame  [Peters  (1999)].  Currently,  to  describe  the  response  of  flame  speed  to  unsteady 
stretch,  expected  in  turbulent  flows,  two  limiting  conditions  are  typically  used.  The  first 
approach  is  to  consider  Markstein  lengths  obtained  in  steadily  stretched  flames  Error!  Reference 
source  not  found,  and  the  second  approach  is  to  consider  vanishing  effect  of  Lewis  number  ( Le ) 
under  overwhelming  presence  of  turbulent  transport  say  in  thickened  flamelet  regimes  Error! 
Reference  source  not  found..  It  is  possible  that  for  wide  conditions  of  interest  in  moderate  to  large 
turbulence  Reynolds  numbers  (Rer  >  100)  and  intermediate  Damkohler  numbers 

(0. 1<  Da  <  10)  the  flame  response  to  stretch  is  somewhat  intermediate.  However,  there  is  little 

systematic  study  on  this  topic,  with  the  exception  of  the  analytical  work  by  Clavin  and  Joulin 
(1997),  although  the  conditions  there  may  not  be  directly  relevant  for  turbulent  flows.  We 
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attempt  to  address  this  problem  by  following/tracking  specific  points  on  the  flame  surface,  by 
displacing  them  with  the  vector  sum  of  their  local  flame  speeds  in  the  surface  normal  direction 
and  local  flow  velocities.  The  stretch  time  history  that  a  flame  particle  carries  can  be  directly 
used  to  answer  fundamental  questions  of  broader  practical  relevance  like  turbulent  flame  speed, 
flame  extinction,  etc.  It  is  noted  that  since  we  only  have  velocity  and  flame  surface  information 
in  a  2D  plane,  we  can  only  track  the  motion  of  a  surface  point  that  is  constrained  to  remain  on 
the  2D  plane.  Adopting  the  definition  and  notation  of  Pope  (1988)  for  surface  points,  the  motion 
of  such  in-plane  flame  particle  can  be  described  by: 

~X,D(0  -  U2D  (x£d  [/],?)+ (x,D  [?p)n£D  (x£d  [(],().  (3,3) 

Here  is  the  position  vector  of  the  flame  particle  on  the  flame  edge  at  time  1.  The  local 
displacement  flame  speed  could  be  modeled  as  Sd  =  S L- 8 M K ,  with  SM  obtained  from  laminar 
expanding  flame  experiments.  Figures  3.3a  and  3.3b  show  the  time  evolution  of  the  flame  edges 
colored  with  local  tangential  strain  rates  Kt  and  stretch  rate  due  to  pure  curvature  K,  respectively 
for  CH4-air  1  atm  2000  rpm  flame  condition.  The  path  tracked  by  a  representative  flame  particle 
is  shown  by  black  circles.  On  the  first  frame  the  flame  particles  are  positioned  on  each  of  the 
flame  edge  points  and  are  allowed  to  evolve  according  to  Eqn.  (3.3).  Since  the  edge  is  located 
only  at  discrete  points  on  the  grid,  it  might  happen  that  the  flame  particles  do  not  land  up  exactly 
on  the  edge  at  the  subsequent  instant.  In  that  case  the  edge  point,  which  is  at  the  shortest  distance 
from  the  new  flame  particle  position,  is  considered  to  be  the  updated  position.  The  Ks  and  K,  time 
histories  from  each  flame  particle  are  subsequently  auto-correlated.  The  integral  of  the  ensemble 
averaged  autocorrelation  function  yields  the  integral  time  scale  for  Kt  and  Kc.  These  time  scales 
are  averaged  over  all  the  runs  over  the  same  conditions  to  yield  the  final  persistence  time  scales 
tp  K  and  tp  K  of  K,  and  K,  respectively.  For  example: 


Kc(^-F2d  [t],t)Kc(X2D[t  +  T],t  +  TyrJ,  (3.4) 

where  the  bar  and  angular  brackets  denote  time  and  ensemble  average  respectively.  These  are 
shown  in  Fig.  3.3c  for  all  experimental  conditions.  We  find  that  for  all  conditions  under  study 
tp  K  remains  rather  invariant  within  a  range  of  400-600  ps.  However,  the  zp  K  shows  a  much 

higher  degree  of  sensitivity  in  response  to  the  change  in  urms,  pressure  and  hence  ReT,  and  varies 
from  about  xPK  -800 jus  for  the  CH4,  1  atm,  4000  rpm  case  to  about  zPK  ~  1 .Sms  for  the 

C4H  10,  5  atm,  2000rpm  case.  In  general  considering  all  cases  under  study,  the  Zp  ^  is  greater  than 

zp  K  by  nearly  a  factor  of  2.  This  suggests  that  strain  rate  due  to  pure  curvature  might  be  a  more 

persistent  quantity  than  that  due  to  tangential  strain  in  the  current  configuration  and  conditions. 
Even  without  flame  particle  tracking,  from  Figs.  3.3a  and  b  we  qualitatively  observe  that  the 
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geometry  of  the  flame  edges  remains  correlated  with  the  evolution  in  time  for  both  cases.  Local 
curvature  being  a  measure  of  the  local  geometry  thus  quantitatively  emerges  as  a  rather  persistent 
quantity  for  the  flames  under  consideration.  Though  these  persistence  time  scales  are  obtained 
from  2D  flame  particle  tracking,  it  is  to  be  recognized  that  experimentally  tracking  the  motion  in 
3D  is  extremely  challenging  if  not  currently  impossible,  given  the  necessity  of  simultaneous 
volumetric  imaging  of  flow  and  reactive  isosurfaces  with  high  time  resolution.  Hence  the 
representative  persistence  time  scale  information  of  tangential  strain  rate  and  curvature  gleaned 
from  the  2D  motion  of  the  surrogate  flame  particle  remains  to  be  of  substantial  interest. 
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Fig.  3.3  :  Evolution  of  flame  edge  with  time,  coloured  by  (a)  K,  (b)  K,  for  CH4-air  0=0.9, 
2000rpm,  p=latm;  (c)  Persistence  time  scales  of  Kc  and  K,  over  ail  experimental  conditions.  CH4- 
air  was  used  in  the  first  five  cases 

4.  Effect  of  Turbulence  on  Ignition  Kernel: 

Ignition  of  a  combustible  mixture  by  a  stimulus  kernel  is  of  relevance  to  many  terrestrial  and 
extraterrestrial  phenomena  and  applications,  such  as  the  various  modes  of  engine  operation 
[Reitz  (2013)],  prevention  of  accidental  explosions  [Kissel  et  al  (1973)],  and  initiation  of  the 
supernova  explosion  [Woosley  et  al  (2011)].  Frequently,  such  ignition  occurs  in  flows  that  are 
highly  turbulent,  involving  a  wide  range  of  time  and  length  scales.  Previous  studies  have 
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advocated  that  turbulence  renders  ignition  more  difficult  [Ballal  (1977),  Kaminiski  et  al  (2000), 
Bradley  et  al  (2004),  Ahmed  et  al  (2006),  Shy  et  al  (2010)],  based  on  the  notion  that  turbulence 
increases  the  dissipation  rate  of  the  deposited  kernel  energy  before  an  embryonic  flame  either  has 
the  time  or  is  aerodynamically  favorable  to  develop. 

Such  an  argument,  however,  does  not  take  into  consideration  the  evolution  and  dynamics  of  the 
structure  of  the  nascent  flame  kernel  after  it  is  formed.  To  appreciate  these  influences,  and  to 
facilitate  exposition  of  the  rationale  of  the  experimental  investigation  to  be  presented  later,  it  is 
first  noted  that  the  burning  rate  and  consequently  the  extinction  propensity  of  a  laminar  flame 
segment  of  surface  area  A,  subjected  to  generalized  aerodynamic  stretching  which  accounts  for 
the  collective  effects  of  flow  nonuniformity,  flame  curvature  and  flame  unsteadiness,  is 
characterized  by  the  stretch  rate  K  =  d(\nA)/dt  [Law  (2006)].  Analysis  in  the  limit  of  large 
activation  energy  Ea  for  an  assumed  one-step  overall  reaction  typically  yields  a  relation  [Ronney 
and  Shivashinsky  (1989),  Sun  and  Law  (2000),  Law  (2006)  ]  having  the  functional  form: 

52ln  s2  =  -  one  (4.1) 

where  s  is  the  local,  stretched  flame  speed  scaled  by  the  steady  unstretched  flame  speed,  k  is  the 
stretch  rate  K  scaled  by  the  laminar  flame  time,  and  a  is  commonly  referred  to  as  the  Markstein 
number,  which  indicates  the  sensitivity  of  the  flame  responses,  such  as  the  propagation  speed,  to 
stretch.  Various  theoretical  expressions  for  a  have  been  derived  [Pelce  and  Clavin  (1982), 
Betchold  and  Matalon  (2001),  Matalon  and  Matkaowsky  (1982),  Clavin  and  Grana-Otero 
(2011)],  accounting  for  effects  of  diffusive  transport  and  thermal  expansion,  the  temperature 
dependence  of  the  transport  properties,  mixture  composition,  general  reaction  orders,  and 
recently  a  modeled  two-step  kinetic  scheme.  These  results  show  that,  depending  on  the  definition 
of  the  flame  location  and  specification  of  the  chemical  kinetics,  the  specific  form  for  ac an  be 
different  for  flame  curvature  and  flow  straining.  The  dominant  influence,  however,  is  the  extent 
that  the  mixture’s  effective  Lewis  number  ( Le )  deviates  from  unity,  as  quantified  by  the 
parameter  (Le  -  1),  where  Le  is  defined  as  the  ratio  of  its  thermal  diffusivity  to  the  controlling 
mass  diffusivity.  This  is  demonstrated  for  example  by  the  expression  [Matalon  et  al  (2003)], 

a  raA(x)  ,  Ze(Le-l)  r^A(x). 

a  = -  - -dx  + - -  —LA.  in 

(cr-1)*'1  x  2(<t-1)  ^  x 


cr-i 


ax 


(4.2) 


where  a  is  the  thermal  expansion  ratio  across  the  flame,  A(x)  the  normalized  thermal 
conductivity  as  a  function  of  the  normalized  temperature  x,  Ze  =  Ea(T\r  Tu)/RTj,~  the  Zel’dovich 
number,  R  the  universal  gas  constant,  and  Th  and  Tu  the  burned  and  unburned  gas  temperatures, 
respectively. 

Two  observations  can  be  made  from  the  above  theoretical  results.  First,  by  examining  the 
linearized  form  of  Eq.  4.1,  s  ~  1  -  aic  it  is  apparent  that  the  trend  of  the  flame  response  is 
qualitatively  affected  by  the  sign  of  a  combined,  diffusivity-affected  stretch  term,  A  =  (Le  -  l)/c, 
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weakened  for  A  >  0,  and  strengthened  otherwise.  Second,  as  a  consequence  of  the  nonlinear 
feedback  between  transport  and  chemical  heat  release  [Sun  and  Law  (2000)],  Eq.  4.1  exhibits  a 
dual-solution  response  for  A  >  0,  indicating  the  potential  of  extinction  at  the  associated  turning 
point. 


Fig.  4.1:  Flame  speed  versus  flame  radius  for  (A)  lean  H2/air  at  (f>-  0.18,  and  (B)  rich  HVair  at  (j>  = 
5.1,  at  different  ignition  voltages  and  turbulent  levels.  The  flame  radius  in  turbulent  cases  is  defined 

as  R  =  yj A/ n  ,  where  A  is  the  area  enclosed  by  the  flame  edge  tracked  from  the  Schlieren  image. 
The  spark  gap  distance  dg  is  0.30  nun  for  A  and  0.58  nun  for  B 


Applying  the  above  concept  to  the  sustenance  of  a  spark-ignited,  continuously-expanding 
spherical  flame  in  quiescence,  it  is  first  recognized  that  since  the  stretch  rate  is  given  by  K  = 
2d(\nR)/dt  >  0,  the  flame  burning  rate  is  governed  by  the  sign  of  (Le  -  1).  It  is  then  obvious  that, 
if  enough  ignition  energy  is  supplied  to  establish  a  flame  kernel,  then  a  Le  <  1  flame  would  be 
strengthened  by  stretch  and  thus  continuously  propagate,  albeit  with  progressively  reduced 
velocity,  until  it  reaches  the  adiabatic  planar  flame  limit  as  R  — »  oo.  However,  for  a  Le  >  1  flame, 
since  the  burning  rate  is  weakened  by  stretch,  even  it  can  be  initially  established,  the  ignition 
energy  may  not  be  sufficient  to  sustain  flame  propagation  when  all  this  energy  is  dissipated 
before  R  has  reached  a  critical  flame  radius,  Rc,  leading  to  extinction.  Thus  if  the  ignition  energy 
is  slightly  larger  than  that  barely  sufficient  to  sustain  flame  propagation  at  Rc,  then  the 
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propagation  velocity  will  exhibit  a  minimum  at  Rc,  and  increase  thereafter  to  attain  the  adiabatic 
planar  flame  limit  as  R  — >oo.  Consequently,  we  expect  the  flame  speed  to  be  small  at  the  state  of 
Rc.  Indeed,  unsteady  analyses  and  simulations  of  the  initiation  of  the  spherical  flame  [Frankel 
and  Shivashinsky  (1984)]  have  shown  that  Rc  corresponds  well  to  the  Zel’dovich’s  flame  ball 
radius,  Rz,  with  a  vanishing  propagation  flame  speed  [Zel’dovich  et  al  (1985)]. 

The  above  distinctively  different  behaviors  are  respectively  shown  in  Figs.  4.1  A  and  4. IB,  for 
the  instantaneous  propagation  velocity  dR(t)/dt,  for  lean  and  rich  FF/air  mixtures  whose  Le's  are 
respectively  <  1  and  >  1;  the  experimental  specifications  are  given  later  when  presenting  the 
results  in  turbulence.  It  is  seen  from  the  solid  lines  that,  for  the  quiescent  situation  of  no 
turbulence,  Mrms  =  0,  an  Le  (=  0.3)  <  1  flame  will  continuously  propagate,  with  monotonically 
decreasing  velocity,  as  long  as  sufficient  ignition  energy,  varied  by  the  spark  discharge  voltage 
f/ign,  is  supplied.  However,  an  Le  (=  2.3)  >  1  flame  will  extinguish  if  the  ignition  energy  is  not 
sufficient  (f/ign  =  80  V),  but  with  sufficient  ignition  energy  (C/ign  =  120  V)  will  continuously 
propagate  after  having  attained  a  minimum,  critical  radius. 

Let  us  now  consider  ignition  in  turbulence,  and  for  the  sake  of  clarity  focus  on  the  critical  state  at 
which  the  ignition  energy  is  just  barely  sufficient  to  initiate  an  expanding,  Le  >  1  flame  in 
quiescence.  It  is  clear  that,  in  the  presence  of  turbulence,  the  otherwise  spherical  flame  with  a 
uniform,  positive,  curvature  will  be  wrinkled,  leading  to  an  embryonic  flame  structure  with  both 
locally  positive  and  negative  curvatures.  Furthermore,  local  flows  with  extensive  and 
compressive  strain  rates,  with  K  >  0  and  <  0,  can  also  be  induced.  The  net  effect  is  that  the  flame 
surface  will  be  locally  subjected  to  additional  A  >  0  and  <  0  stretch  effects.  The  key  point  to 
note  here  is  that  since  the  flame  in  quiescence  is  already  at  the  incipient  state  of  extinction  over 
its  entire  surface,  any  flame  segments  that  are  subjected  to  additional  A  >  0  effect  in  turbulence 
are  still  susceptible  to  extinction.  However,  the  burning  intensity  of  the  flame  segments  that  are 
subjected  to  A  <  0  effect  will  be  increased,  moving  them  away  from  incipient  extinction  and 
consequently  can  collectively  serve  as  local  ignition  sources  to  sustain  the  global  flame  structure. 

With  the  above  expository  anticipation,  the  present  work  then  aims  to  explore  such  a  possibility, 
which  may  alter  the  traditional  view  on  the  criteria  of  ignition  in  turbulent  flows,  and  as  such 
highlight  the  necessity  to  incorporate  the  dynamics  of  the  embryonic  flamelet  structure  into  the 
description  and  prediction  of  ignition. 

We  undertook  a  well-controlled  experimental  approach  to  categorically  study  the  effects  of 
turbulence  on  an  ignition  kernel.  Tests  were  conducted  in  the  constant-pressure,  optically 
accessible,  vessel,  detailed  in  the  section  of  experimental  details.  The  ignition  system  was  similar 
to  those  used  in  automotive  engines,  with  a  slight  modification  by  replacing  the  spark  plug  by 
two  tungsten  wires  of  250-pm  diameter  and  centered  axially.  The  spark  is  generated  by 
discharging  a  33  pF  capacitor  through  an  ignition  coil  with  a  134:1  turn-ratio.  Voltage  (C/ign) 
across  the  capacitor  was  made  variable  from  50  V  to  170  V,  allowing  the  ignition  energy  to  vary. 
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In  designing  the  matrix  of  the  experimental  investigation,  it  was  noted  that  previous  works  on  the 
interaction  of  turbulence  and  ignition  kernel  were  limited  to  inhibited  ignition,  ostensibly  due  to 
the  small  deviations  in  Le  from  unity.  Therefore  we  extended  the  experiments  to  cover  a  wide 
range  of  Le ,  especially  focusing  on  mixtures  with  Le  sufficiently  greater  than  unity  such  that 
ignition  could  be  limited  by  the  extent  of  stretch  that  the  flame  kernel  experiences.  At  the  same 
time  we  also  considered  mixtures  with  Le  ~  1  or  Le  <  1  to  provide  a  complete  description  of  the 
phenomena.  We  selected  H2/O2  mixtures  for  most  of  the  investigation  because  it  allows  flexible 
variation  of  Le  due  to  the  drastic  difference  between  the  molecular  weights  of  hydrogen  and 
oxygen.  This  can  be  achieved  by  varying  the  H2/O2  ratio,  namely  the  equivalence  ratio  (f)  of  the 
mixture  (defined  to  be  the  fuel/oxygen  ratio  relative  to  the  stoichiometric  mixture),  thereby 
changing  the  diffusivity  of  the  deficient,  hence  controlling,  reactant.  Furthermore,  since  the 
H2/O2  oxidative  chemistry  is  relatively  simple  and  also  reasonably  well  established,  the  present 
data  can  be  usefully  adopted  in  further  computational  studies  on  issues  of  turbulence-chemistry 
interaction  in  general  and  of  this  practically  relevant  problem  in  particular. 

We  first  demonstrate  the  possible  enhancement  of  ignition  by  turbulence  through  direct,  time- 
resolved  Schlieren  images,  using  a  fuel  rich  FF/air  mixture  of  (f>=  5.1,  which  has  a  large  Le  of  2.3 
since  the  controlling,  deficient  reactant  is  Oo,  while  the  thermal  diffusivity  is  controlled  by  that 
of  H2.  Figs.  4.2A  and  4.2B  respectively  show  a  successful  and  a  failed  ignition  event  in 
quiescence  («rrns  =  0),  with  an  ignition  voltage,  f/]gn  =  120  V  and  a  smaller  value  of  80  V.  This 
result  therefore  agrees  with  the  earlier  discussion  on  the  dependence  of  ignition  on  the  energy 
input.  The  measured  velocities  of  these  images,  together  with  those  for  the  lean  mixture  of  </>  = 
0.18,  with  Le  =  0.3,  are  those  shown  in  Figs.  4.  IB  and  4.1A,  respectively. 

If  we  next  maintain  the  ignition  voltage  constant  at  80  V  but  progressively  increase  the 
turbulence  level,  Figs.  2B-E  show  that  the  structure  of  the  flame  kernel  is  changed  from  a 
positively  stretched  spherical  surface  to  a  multitude  of  wrinkled  flamelets  of  both  positive  and 
negative  curvatures,  and  that  ignition  is  achieved  at  u „  =  2.9  m/s.  This  result  therefore  supports 
the  notion  that  for  Le  >  1,  a  combustible  mixture  can  be  ignited  in  turbulence  with  an  ignition 
energy  that  is  not  enough  to  ignite  the  same  mixture  in  quiescence. 

We  have  extensively  mapped  out  the  ignition  boundary  for  the  (/>  =  5.1  mixture  in  terms  of  the 
turbulence  intensity  and  the  ignition  energy,  summarized  in  Fig.  3A.  It  is  found  that  the 
minimum  discharge  voltage  required  for  successful  ignition  in  the  presence  of  moderate 
turbulence  («rms  =  2.9  m/s)  is  65  V  (indicated  by  blue  circles  filled  with  green  color  in  figure), 
which  is  substantially  reduced  from  the  voltage  required  for  ignition  in  quiescence  (100  V). 
Based  on  the  relation  Elgn  =  Cf/Ign2,  where  Eign  and  C  are  the  discharge  energy  and  capacitance, 
respectively,  this  result  implies  that  the  minimum  ignition  energy  can  be  lowered  by  a  few 
factors  with  the  presence  of  turbulence. 
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Fig.  4.2:  Sequential  Schlieren  images  of  flame  kernel  development  for  H2/air  at  <f>  =  5.1  (Le  ~  2.3), 
at  different  ignition  voltages  and  turbulent  levels.  Schlieren  imaging  signal  is  proportional  to  the 
density  gradient  of  the  flow  field,  a  good  indicator  of  the  flame  front.  The  view  for  each  image  is  65 
mm  x  65  mm.  t  is  time  after  discharge. 


Having  demonstrated  the  facilitating  effect  of  turbulence  on  the  Le  >  1,  rich  Hi/air  mixtures,  it  is 
necessary  to  investigate  the  response  of  lean  mixtures,  whose  Le  is  less  than  unity  because  the 
controlling,  deficient  reactant  is  now  H2  and  the  thermal  diffusivity  is  dominated  by  those  of  02 
and  N2.  As  anticipated  earlier,  since  a  flame  kernel  with  Le  <  1  is  always  facilitated  by  positive 
stretch  in  quiescence,  ignition  is  expected  to  be  successful  as  long  as  a  flame  kernel  can  be 
established.  Turbulence,  in  this  case,  increases  the  dissipation  rate  of  the  deposited  kernel  energy 
and  as  such  renders  ignition  progressively  more  difficult  with  increasing  turbulence  intensity, 
eventually  leading  to  extinction.  Our  results  on  lean  H2/air  mixtures  (</>=  0.12  -  0.2,  Figs.  3B  and 
3C)  confirm  this  mechanistic  interpretation  and  are  also  consistent  with  the  conclusions  of 
previous  studies  which,  as  noted  earlier,  were  mostly  based  on  Le  ~  1  mixtures. 
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Fig.  4.3:  Ignition  test  results  plotted  against  turbulence  intensity  (unm )  and  discharge  voltage  (Ujgn. 
)  or  mixture  reactivity,  represented  by  the  H2/O2  ratio  or  the  amount  of  inert  gas.  ( A )  H2/air  at 
equivalence  ratio  ( (f>)  =  5.1  ( B )  H2/air  at  (f>  =  0.12  (C)  H 2/air  for  fixed  Ujgn  =  160  V.  (D)  H2/02/He 
at  (f>  =  1.0,  02/(02+He)  Vol.%  =  8.0%  (E)  H2/02/C02  at  0=  1.0,  02/(02+  C02)  Vol.%  =  10.3%  (F) 
H2/O2/CO2  at  f  -  1.0  for  fixed  U,gn  =  160  V  (G)  H2/02/Ar  at  <f>  =  1.0,  02/( 02+  Ar)  Vol.%  -  5.0% 
(H)  H2/02/Ar  at  <j)-  1 .0  for  fixed  Uign  =  160  V.  (I)  n-C4H  1(/air  at  (ft  =  0.7  (J)  n-C4H  1(/air  at  (f>  =  2.2 
(K)  n-C4H  10/ air  for  fixed  Ulgn  -  160  V.  Initial  pressure  and  temperature  for  all  tests  are  1  atm  and 
298  K.  The  Le  shown  in  figure  is  defined  to  be  the  effective  Lewis  number  of  a  combustion  mixture, 
based  on  theory  in  [24].  The  spark  gap  distance  dg  is  0.58  mm  for  case  A,  0.30  mm  for  cases  B-H, 
and  0.80  mm  for  cases  I-K.  Here  we  identify  three  stages.  Failed  Ignition  (red  crosses ):  where 
ignition  kernel  fails  to  grow  to  a  propagating  flame,  either  due  to  lack  of  ignition  energy  in 
quiescent  environment  or  faster  dissipation  by  turbulence.  Successful  Ignition  (green  dots):  where 
ignition  kernel  containing  sufficient  energy  grows  to  a  flame  with  or  without  the  presence  of 
turbulence.  Ignition  Facilitated  by  Turbulence  (blue  open  circles  filled  with  green  color):  where 
ignition  kernel  grows  to  aflame  in  the  presence  of  turbulence,  but  fails  to  do  so  in  quiescence  if 
supplied  with  the  same  ignition  energy. 


We  also  need  to  rule  out  the  possibility  that  the  controlling  ignition  chemistry  could  be  different 
for  fuel  rich  and  fuel  lean  mixtures,  causing  the  observed  phenomena.  To  remove  such  a 
possibility,  we  manipulated  the  values  of  Le  by  changing  the  inert  bath  gas  from  No  to  (He,  Ar, 
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CO2),  with  the  fuel/oxygen  ratio  fixed  such  that  the  controlling  lean-vs.-rich  chemistry  is  not 
affected.  It  is  seen  (Figs.  3D-H)  that  the  facilitating  effect  of  turbulence  is  again  manifested  only 
for  Le  >  1  mixtures,  and  that  the  effect  can  be  flipped  from  inhibiting  to  facilitating  solely  by 
changing  the  inert  bath  gas  while  (f>  is  fixed  at  unity.  In  particular,  turbulence  facilitates  ignition 
with  He,  a  light  inert  which  substantially  increases  the  thermal  conductivity  and  thus  Le\  while 
with  heavier  inerts  such  as  CO2  and  Ar,  Le  becomes  less  than  or  near  unity  and  as  a  result 
ignition  is  inhibited.  These  extensive  sets  of  experiments  therefore  rule  out  the  possibility  that 
the  distinct  turbulence  effects  for  fuel  lean  and  fuel  rich  conditions  are  due  to  different  chemical 
kinetics. 

Finally,  in  order  to  provide  even  further  substantiation  of  the  phenomena  and  concept  advanced 
herein,  and  for  fuels  other  than  H2,  we  have  flipped  the  lean-vs.-rich  aspect  of  the  mixture  by 
using  77-butane  (77-C4H10)  as  the  fuel,  whose  Le  are  greater  and  smaller  than  unity  for  lean  and 
rich  mixtures  because  of  the  substantially  larger  molecular  weight  of  77-butane  relative  to  that  of 
oxygen.  Thus  lean  (rich)  72-butane/air  mixtures  should  exhibit  behavior  similar  to  those  of  rich 
(lean)  H2/air  mixtures,  respectively. 

Figs.  4.3I-K  summarizes  the  ignition  test  results  for  a  lean  and  a  rich  77-butane/air  mixture,  whose 
Le’ s  are  2.1  and  0.9  respectively.  It  is  seen  that  ignition  is  indeed  facilitated  for  the  Le  >  1,  lean 
mixture.  The  above  extensive  results  therefore  demonstrate  that  ignition  enhancement  by 
turbulence  for  mixtures  with  large  Le  is  of  a  general  nature,  irrespective  of  fuel,  equivalence  ratio 
and  inert. 

The  above  result  unambiguously  shows  the  role  of  Le  in  propagation  of  an  embryonic  ignition 
kernel  in  a  turbulent  environment.  More  categorically,  it  shows  that  for  Le  >  1  conditions, 
turbulence  enhances  the  growth  of  the  kernel  in  presence  of  turbulence  as  opposed  to  that  in  the 
case  of  Le  <  1,  where  turbulent  mixing  or  dissipation  simply  blows  away  the  kernel  before  it 
develops  into  a  flame.  This  immediate  raises  a  question  that  if  dissipation  is  not  important  at  all 
for  the  Le  >  1  conditions.  To  answer  the  question  we  studied  such  conditions  of  Le  >  1, 
primarily  with  lean  72-butane  (72-C4H10)  mixtures,  with  increasing  turbulent  intensities.  We  found 
for  the  ignition  voltage/energy  where  turbulence  assisted  ignition  is  possible,  i.e.,  even  though 
the  ignition  fails  at  the  laminar  condition  with  certain  level  of  turbulence  the  kernel  grows  into  a 
flame,  there  exists  an  upper  boundary  of  turbulence  intensity,  beyond  which  turbulent  suppresses 
ignition.  These  results  and  the  complete  map  for  lean  72-butane  (72-C4H10)  are  shown  in  Fig.  4.4, 
which  shows  the  non-monotonic  response  of  the  ignition  kernel  with  the  turbulence  intensities. 
This  is  caused  by  the  increased  dissipation  with  the  higher  turbulence  intensity,  which  overcomes 
the  benefits  of  locally  negative  stretched  flamelets  whose  strong  propagation  for  Le  >  1 
mixtures  caused  the  turbulence  assisted  ignition  possible  for  lower  levels  of  turbulence  intensity. 
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Fig.  4.4:  Ignition  map  for  n-C4H I0/air  at  (f>  =  0.7.  Here  we  identify  the  ignition  window  with  the 
presence  of  three  stages.  Failed  Ignition  (red  crosses):  where  ignition  kernel  fails  to  grow  to  a 
propagating  flame,  either  due  to  lack  of  ignition  energy  in  quiescent  environment  or  faster 
dissipation  by  turbulence.  Successful  Ignition  (green  dots):  where  ignition  kernel  containing 
sufficient  energy  grows  to  a  flame  with  or  without  the  presence  of  turbulence.  Ignition  Facilitated 
by  Turbulence  (blue  open  circles  filled  with  green  color):  where  ignition  kernel  grows  to  aflame  in 
the  presence  of  turbulence,  but  fails  to  do  so  in  quiescence  if  supplied  with  the  same  ignition 

In  summary,  while  turbulence  is  usually  believed  to  suppress  ignition  due  to  the  enhanced 
dissipation  of  localized  ignition  energy,  we  have  experimentally  demonstrated  that  it  can  actually 
facilitate  ignition  under  conditions  in  which  ignition  is  limited  by  the  difficulty  of  the  flame 
kernel  to  transition  into  an  expanding  flame.  This  is  possible  through  an  embryonic  flamelet 
structure  consisting  of  segments  subjected  to  both  positive  and  negative  stretch,  while  quiescent 
ignition  generates  only  positive  stretch  through  the  positive  curvature  over  the  entire  expanding, 
smooth  flame  surface.  Such  an  understanding  is  of  both  practical  and  scientific  significance.  For 
example,  explosion  tests  in  quiescence  may  underestimate  the  risk  of  accident;  while  engine 
flows  can  also  be  optimized  to  reduce  misfire  for  ultra-lean,  clean,  fuel-efficient  operations, 
recognizing  that  engine  fuels  are  large  hydrocarbons  such  that  their  lean  burning  corresponds  to 
Le  >  1  situations.  Furthermore,  while  studies  of  the  supernova  explosion  have  always  assumed 
the  initial  existence  of  a  flame  that  subsequently  transitions  to  a  detonation  wave,  perhaps  it 
would  be  of  interest  to  also  investigate  situations  for  which  the  establishment  and  sustenance  of 
such  a  flame  may  not  be  possible  in  the  first  place,  recognizing  the  fact  that  the  Le  for 
supernovae  is  exceedingly  large  [Glazyrin  et  al ],  of  the  order  of  104. 
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